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Abstract 

We consider a single level quantum dot interacting with a phonon 
mode and weakly coupled to metallic leads which are subjected to a time- 
dependent gate voltage. Electrical conductance, thermopower, and figure 
of merit are investigated in detail using a Tien-Gordon-type rate equation. 
In the presence of the microwave field, the electrical conductance exhibits 
extra peaks whose height is controlled by the magnitude of the microwave 
field. Furthermore, the oscillation of the thermopower increases in the 
presence of the time-dependent gate voltage or the electron-phonon in- 
teraction. Influence of the electron-phonon interaction, microwave field, 
temperature, and Coulomb interaction on the figure of merit is also stud- 
ied. 

1 Introduction 

The study of the thermoelectrical properties of nanoscale devices has attracted 
a lot of attention in recent years theoretically and experimentally [TJ [31 [31 HI |H1 
H M E3 ED 02] ■ The violation of the Wiedemann Franz law [Ml [14] , that the 
ratio of the electrical conductance to the thermal conductance is proportional 
to the temperature, in nanostructures has caused that the nanostructures have 
been considered as good candidates for fabricating thermoelectric devices. The 
thermoelectric efficiency is measured by a dimensionless quantity as figure of 
merit, ZT = S 2 G e T / (n e + n p h) where S stands for the thermopower and G e 
denotes the electrical conductance. K e , and n p h are the electrical and phononic 
thermal conductances, respectively. The Coulomb blockade effect [131 Hj>], dis- 
creteness of the energy levels [TC] , interference effects [TU] and so on are the main 
reasons for high figure of merit in nanoscale devices. The quantum dot (QD) 
as a promising thermoelectric material has been widely studied during recent 
years [3l HI [Til [13l [T5l [16l HZ1 [18l [19l [20l [21] [22l [23l [24l [25]. 

The influence of the time-dependent gate voltage on the thermoelectric prop- 
erties and on the heat transport through QD-based devices is a very interesting 
issue which has been considered very recently. It was predicted that the heat 
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can be conducted by phonons in the low temperature [55] . In addition, it was re- 
ported that the thermopower can be enhanced in the presence of a step- like pulse 
bias [23 [55] . Pei and co-workers [33] studied the transient heat generation in a 
QD by applying a step-like pulse bias. They found that the time-dependent bias 
results in the periodic oscillation of the heat generation. Chi and co-workers [50] 
studied the thermoelectric properties of a QD in the presence of a microwave 
field. They reported that the applying microwave leads to the appearance of 
the extra peaks in the electrical conductance spectrum. 

The coupling between the electronic degree of freedom and vibrational de- 
gree of freedom or the electron-phonon interaction (EPI) is an important issue 
in the transport through nanoscale devices resulting in the novel and interest- 
ing phenomena. The influence of the EPI on the transport properties of the 
molecular transistors and the carbon nanotube quantum dots has been exten- 
sively studied using both the Keldysh nonequilibrium Green function formalism 
and the rate equation approach [SH E31 EH [SHI EH EH IMl HOI SH H2]- 
In addition, the thermoelectric properties of the molecular transistors has been 
analyzed in recent years [TU [331 HH US HB]. However, although the transport 
properties of such devices has been studied in the presence of the EPI and 
the time-dependent voltages [37[ H5J , the influence of the time-dependent gate 
voltage on the thermoelectric properties of the molecular devices has not been 
addressed so far. In this article, we consider a molecule attached to the metal- 
lic leads which are subjected to a microwave field in the non-adiabatic regime. 
We follow the method introduced by Dong and co-workers [48] to obtain the 
formal expressions of the thermopower and electrical conductance. The analyti- 
cal relations for the thermopower and electrical conductance are obtained using 
a Ticn-Gordon-type rate equation. The influence of the EPI, microwave field, 
temperature and Coulomb repulsion on the figure of merit is also examined. The 
model and formalism are presented in the next section. Section 3 is devoted to 
the numerical results, and in the end, some sentences are given as a summary. 

2 Model and formalism 

We consider a single level QD interacting with a dispersionless optical phonon 
mode localized in the QD. The QD is weakly coupled to the normal metal elec- 
trodes which are subjected to a time-dependent gate voltage. The Hamiltonian 
describing the whole system is given as 

H = ^ £ aka{t)c* aka C aklJ + ^ e ° n ° + Un fH + watft + (!) 
aka o 

A^[a f + a]n c + ^[K^c^ck + H.c], 

(7 aka 

where c a ka(c akiJ ) destroys (creates) an electron with wave vector k, spin er in 
lead a. £ a fc CT (£) = e° aka +u a cos^lt denotes the energy levels of the leads composed 
of the rigid shift of the Fermi energies of the leads, £° aka , and the time-dependent 
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energies, it£,(R) = ±V ac induced by a microwave field with the frequency and 
the magnitude V ac . d cr (d\,) is the annihilation (creation) operator in the QD and 
n a = d^da is the occupation operator. U denotes the on-site Coulomb repulsion 
whereas, ui is the phonon energy. A and V a ko describe the electron-phonon 
coupling strength and tunneling between the QD and the leads, respectively. 

The EPI can be eliminated by using a polaronic transformation [49] . H = 
e s He~ s where S = exp((a^ — a) J2 a n <0- The transformation results in decou- 
pling the electronic and phononic subsystems and the renormalization of the 
QD energy levels and the Coulomb repulsion. The transformed Hamiltonian is 
given as 

H = ^2 £akaC f aklT c a kcr + e a n a + Un^ni (2) 

cxka <J 

LOO 1 a + y^JVotkcrX&^da + H.C.] 
aka 

where e„ — e a — A 2 /o>, U = U — 2X 2 /uj, and X = exp{—X/u){a i — a)). In 
the following, we use the wide band approximation so that the tunneling rate, 
= J2 k \ Vakal 2 , is energy-independent. Because the QD-lead coupling is 
assumed to be so weak, « T, where T denotes the temperature, and 
on the other hand, the frequency of the microwave field is so stronger than 
the tunneling rate, a Tien- Gordon- type rate equation [351 ED] can be used to 
describe the behavior of the system. With respect to Eq. , the isolated QD 
can be in the empty state, one-electron state, or doubly occupied state with n 
phonons. These states are shown by PjJ where N and n stand for the number of 
the electrons and phonons, respectively, and obey the following rate equations: 



pn 

±=Y, t r -o P ™ - r X^o] (3a) 

am a 

pn 

a \ ^ r-pmn r>rn , j^mn r>m (j^nm . prim \ nni /qt\ 

"37" — "+" 1 a2<r*2 ~ I 1 aaO ' 1 olg2)^g\ \ 60 ) 
am 

= E[ r -2^ m -rSM (3c) 



pi 



dt 

ama 

where T"^ M describes the transition from N-electron state to M-electron state, 
while the number of phonons is changed from iitom. The transition rates are 
defined as @S] 

r^M = r Q e^ 2 . 9 2 i" l -"i^(4 m -"i( 5 2 )) 2 £ J 3 (u a /nf (4) 

j=—oa 

[f a {e M -e N + (m- n)uj - j9)Q{M - N)+ 
(1 - f a (e N - e M + (n - m)u - jSl))Q{N - M)\ 

where g = A/cj, p = min[m,n], q — max[n,m], and Q(x) is the Heaviside step 
function. N and M denote the number of the electrons in the QD expressed as 
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N = 0, empty state, N — 1, s~n — e a (single electron state), and N = 2, ejy = 
if + i± + U (two electron state). L^(x) is the generalized Laguerre polynomial 
and Jj(x) is the Bessel function of order j. f a (x) = (1 + exp((x — /^ Q )/fcT)) _1 
is the Fermi distribution function of lead a and fi a stands for the chemical 
potential of the lead. Note that the leads are assumed to be the normal metal, 
therefore, the tunneling rate, T Q , is spin-independent. 

Solving Eqs. ([3]) in the steady state (dPfi/dt = 0), the charge and energy 
currents are expressed as 

1 =-^2^ i T aOa P + ( r a CT 2 ~ T aao) P a ~ T a 2a P 2\ (5a) 
anra 

n a _ 1 \rQ nm P™ J- (yQ nm r Qnm\ p n r Qnm p i /ci.'i 

y T l aOa V aa2 ckctO ) (t q2ct 2\ (5b) 



an in 



where r^JJ- 7 ^ describes the energy transported in the tunneling process and 
expressed as 

| oo 

r& = r«e-*yi— »l|(Ll-«l ( ^ ))2 £ J^/fif (6) 

j=-oo 

[(sm - ew + (m- n)uj - jVL)f a {e M -e N + (m- n)ui - jti)@(M - N)+ 
(ejv - e M + (n - m)u - jft)(l - /„(% - e A/ + (n - m)w - jft))6(iV - A/)] 

In the linear response regime, the charge and heat currents are expanded 
as 05] 

J Q = G e A^ + G T AT (7a) 

Q a = MAV + ifAT (7b) 

where G e and Gt are the electrical conductance and thermal coefficient, re- 
spectively. The thermopower is computed as S = —Ay/ AT in the limit of zero 
current. In order to compute the above coefficients, we assume Tl = T + AT, 
Hl = Ef + eAV, = Ef, and Tr = T where Ef denotes the Fermi energy 
of the leads. By expanding the Fermi function of the left lead according to 
f L (e) = /(e) - eAV f '(e) - $f(e - Ef)f(e) where /'(e) = df(e)/d(e) and 
setting I L = 1/2(I L - I R ) and Q L = l/2(Q L - Q R ), one can easily obtain the 
thermoelectric coefficients. For simulation purpose, we set oj = 1 as the energy 
unit. The reported value of the phonon energy is varied from 10[ieV to \QmeV 
in the literatures. In addition, we set n p h = 3kq where kq = ir 2 k 2 /3hT is the 
quantum of thermal conductance [51] and, assume that the energy level of the 
QD is degenerate. 

3 Numerical results and discussions 

Figs, la and lb describe the electrical conductance as a function of the energy 
level of the QD and temperature, respectively. As it is expected, the main peaks 
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of the electrical conductance are located in the resonance energies, i.e., i = 0, 
and e = —U. First, we analyze the behavior of G e in the absence of the EPI. 
By applying microwave field the electrical conductance shows extra peaks at 
the energies e = ±nf2, and s = — U ± nf2 where n is an integer. These side 
peaks are induced by the photon-assisted tunneling channels. Such result has 
been recently reported in Ref. [30] ■ By solving Eqs. (4,7), we have obtained the 
following analytical relation for the electrical conductance 

_ e 2 I* /kT ~ Jn{u L /V) 2 ( ) 

e H2kTZ [y ' £-? l + coah((Ef + nl[l-e)/kT) K > 

e/kT , -{2e+U)/kT\ Jn(uL/ty 2 i 



{e -e/kT + e -(2e + U)/k T) £ _ 



n— — oo 



+ cosh((Ef + nn-s- U)/kTy 



where Z = 1 + 2exp(—e/kT) + exp(—(2e + U)/kT) is the partition function and 
J n {x) is the Bessel function of order n. It is straightforward to show that the 
electrical conductance obeys the following equation for the case V ac = 

e 2 T L 1 e U/kT _|_ e -e/kT 

Ge = J kTZ [ 1 + e £ / fcT + (1 + e {£+u)/kT)2 1 ( 9 ) 

that it is in agreement with the results presented in Ref. 43 . According to 
Eq. ([9]) , the height of the main peaks of the electrical conductance is equal to 
\wr if ml = 0. In the presence of the time-dependent voltage, the height of 
the main peaks decreases proportional to Jq(ul/£1) 2 . Therefore, the ratio of 
the main peak's height in the case ul ^ to ul — can contain the important 
information about the influence of the time-dependent gate voltage on the en- 
ergy levels of the electrodes. Moreover, the appearance of the extra peaks in the 
electrical conductance is explained by Eq. ([5]). The height of these side peaks 
is controlled by J„(iti/fi) 2 so that their height increases with increase of Ul- 
The dependence of the electrical conductance on the temperature is strongly 
dependent on the energy level of the QD, so that the electrical conductance 
is significantly reduced in the resonance energies with increase of temperature. 
The behavior of the electrical conductance in the electron- hole symmetry point, 
e = —U /2, is very interesting. In low temperatures, kT << U , and in the ab- 
sence of the microwave fields, the electrical conductance is zero because the elec- 
trons and holes participate in the current with the same weight. The magnitude 
of the electrical conductance significantly increases in the presence of the time- 
dependent voltage. By analyzing Eq. (j8J) in the low temperature limit, it is found 
that the electrical conductance is equal to % [J_i(u_l/S!) 2 + ^(ul/D,) 2 } so 
the magnitude of G e is increased with increase of Ul . The influence of the pho- 
tons is suppressed in the higher temperatures so that the electrical conductance 
becomes independent of the photon-assisted channels. In the high temperature, 
the various electronic states participate in the charge transport through the QD 
thus, the extra channels induced by photons do not have significant role in the 
charge transport. Now, we study the effect of the EPI on the electrical conduc- 
tance. As one expects, the position of the main peaks is now shifted because 
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of the polaronic shift. It is important to note that the distance between the 
peaks is controlled by the Coulomb interaction considered to be constant in the 
article, U — 2. The height of the peaks is exponentially reduced in the pres- 
ence of the EPI. A careful analysis reveals that the reduction is proportional 
to the EPI strength as e~ 9 . Such decrease was also observed in the Kondo 
regime [52] owing to the existence of factor e~ 9 in the density of the states of 
the QD. Here, the reduction is caused due to the exponential suppression of the 
tunneling between the QD and the leads in the presence of the EPI. Results 
also show that the EPI results in the reduction of the electrical conductance, 
more specifically, in high temperatures. However, the photon-phonon-assisted 
electrical conductance is more than the elastic case in low temperatures. The 
decrease of the phonon-assisted electrical conductance comes from the fact that 
the electron configurations with few phonons are now partially populated but 
their participations in the conductance are suppressed by |L^ >0 (<7 )| . 

Fig. 2 describes the thermopower as a function of the energy level of the 
QD. The sawtooth oscillation of the thermopower is risen due to the change of 
the electron number in the QD. By solving Eqs. (4,7), the following analytical 
relation for the thermopower is obtained in the case A = 



When ul — 0, the thermopower has three zeros located in the resonance ener- 
gies and electron- hole symmetry point. In resonance energies the temperature 
gradient cannot produce the net current and, as a result, the thermopower be- 
comes zero. In the electron-hole symmetry point, both electrons and holes carry 
the charge and energy with the same weight. Although they carry the energy 
in the same direction, they carry the charge in the opposite directions so that 
no net current is produced. The sign of the thermopower shows the kind of 
charge participated in the transport so that S > is for holes, while S < 
is for electrons. In the energies more than the resonance energies electrons in- 
jected from the hotter lead (the left lead) create the current, whereas for the 
right hand side of the electron-hole symmetry point, holes create the current. 
Such behavior was previously reported for single QD and double quantum dot 
systems [TBJ [531 HI] • I n the presence of the microwave fields and EPI, the os- 
cillations of the thermopower increase because of the phonon-photon-assisted 
tunneling channels. However, the magnitude of the thermopower decreases pro- 
portional to J„(ul/51) 2 . Moreover, oscillations become faster with increase of 
ul, because the energy levels of the leads oscillate faster. The slope of the ther- 
mopower is roughly estimated like 1/T which is in agreement with the results 
presented in Ref. [43] for single molecule devices in the sequential tunneling 
regime. In the energies near the chemical potential of the leads, the sign of 
the thermopower changes more because the phonons and photons have more 
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meaningful role in the charge transport from the leads to the QD. 

The dependence of the figure of merit on the energy level of the QD and 
temperature is shown in figs. 3a and 3b, respectively. The figure of merit 
is zero in the resonance energies and electron-hole symmetry points because of 
5 = 0. The EPI and microwave fields result in the reduction of the magnitude of 
the figure of merit because they modulate the tunneling process between the QD 
and the electrodes so that the thermopower and electrical conductance decrease. 
However, the number of peaks of ZT significantly increases in the presence 
of the microwave fields or the EPI owing to the photon- or phonon-assisted 
channels. Indeed, the time-dependent gate voltage results in the increase of the 
thermoelectric efficiency of the device in energies in which the device cannot 
work as a thermoelectric device without the microwave fields. The same result 
was reported by Chi and co-workers [30]. In the low temperatures, the figure of 
merit has a maximum (fig. 3b) because the Coulomb interaction controls the 
transport process. With increase of temperature resulting in the reduction of 
the Coulomb interaction effect, the figure of merit is reduced. The position of 
the peak is slightly shifted in the presence of the EPI or the microwave fields. 

Fig. 4a shows the dependence of ZT on the strength of the EPI and the 
magnitude of the microwave field. Results show that applying the microwave can 
partly compensate the negative effect of the EPI on the figure of merit. It comes 
from the photon-assisted channels resulting in the increase of the thermopower 
and electrical conductance even in the presence of the EPI. Therefore, the time- 
dependent gate voltage can be used to improve the thermoelectric efficiency of 
the molecular devices in the presence of the EPI. The influence of the Coulomb 
interaction on the figure of merit is analyzed in fig. 4b. In the range of U = — 2e, 
(e = —0.5), ZT is zero because of the electron-hole symmetry point. Indeed, 
ZT approaches zero when the electrons and holes participate in the charge and 
energy transport. In the both weak and strong Us, ZT takes the finite values. In 
these cases, electrons carry the charge and energy through the system because 
the QD is in the two-electron state or one-electron state, respectively. It is 
interesting to note that ZT takes the reasonable values in high As only in strong 
Us. It comes from the fact that the empty state ore two-electron state will be 
populated in strong A and nearly weak U, so that holes can also participate 
in the transport resulting the decrease of the thermopower. Therefore, ZT 
becomes very small in weak U and strong A. 

4 Summary 

In this article, we have investigated the thermoelectric properties of the molec- 
ular devices in the presence of the time-dependent voltage by means of a Tien- 
Gordon- type rate equation formalism. Analytical expressions of the thermopower 
and electrical conductance are obtained. Results show that the microwave field 
results in the appearance of the extra peaks in the conductance spectrum whose 
height is controlled by the magnitude of the microwave. Electron-phonon inter- 
action results in the reduction of the figure of merit. However, the microwave 
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field can slightly compensate the negative effect of the EPI. We also find that 
the figure of merit takes reasonable values in the strong Coulomb interaction if 
the EPI is very strong. 
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Figure captions 

Figure 1: Electrical conductance (a) as a function of the energy level of the 
QD, (b) as a function of temperature. Parameters are f2 = 0.5, T = O.Olw, and 
U = 2lo. kT = O.lin fig. la and e = -w in fig. lb. 

Figure 2: Thermopower versus the energy level of the QD. Parameters are 
the same as fig. 1 except e = — 0.5w. 

Figure 3: Dependence of ZT on (a) energy level of the QD, (b) temperature. 
Parameters are the same as fig. 1. e = — 0.5uj in fig. 3b. 

Figure 4: Color map of ZT. Parameters are the same as fig. 3. V ac = uj in 
fig. 4b. 
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Figure.l 
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